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After deriving a general correspondence between linear response correlation functions in graphene 
with and without applied uniaxial strain, we study the dependence on the strain modulus and 
direction of selected electronic properties, such as the plasmon dispersion relation, the optical con- 
ductivity, as well as the magnetic and electric susceptibilities. Specifically, we find that the dispersion 
of the recently predicted transverse plasmon mode exhibits an anisotropic deviation from linearity, 
thus facilitating its experimental detection in strained graphene samples. 
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I. INTRODUCTION 

Graphene is a truly two-dimensional (2D) electronic 
system, based on an atomically thin carbon honey- 
comb lattice^. Low-energy quasiparticles can be de- 
scribed as massless Dirac fermions, with a cone disper- 
sion relation in reciprocal space around the so-called 
Dirac points K, K', and a linearly vanishing density 
of states (DOS) at the Fermi leveP''^. Such a linear 
spectrum and reduced dimensionality yield remarkable 
behaviors already in the non-interacting limit of sev- 
eral electronic properties of graphene. These include, 
inter alia, the reflectivitj*^, the optical conductivitj«^""— , 
the plasmon dispersion relationjiS^J^, as well as a newly 
predicted transverse electromagnetic mod©i^, which is 
characteristic of a 2D system with a double band struc- 
ture, such as graphene. Moreover, the relevance of 
lossless plasmons in graphene in the infrared frequency 
range has been emphasized, with possible applications in 
nanophotonics^^'^^. These properties can be extracted 
from the study of the appropriate correlation functions 
within linear response theoryii"— . 

Within the Dirac approximation, rotational invariance 
implies thet the current-current correlation function may 
be decomposed in a longitudinal and a transverse contri- 
bution, the former being related to the density-density 
correlation function via the continuity equatio n^^i^" . 
Moreover, in the case of massless Dirac fermions, it has 
been observed that the current-current response is simply 
proportional to its pseudospin-pseudospin counterpart^". 
On the other hand, the magnetic susceptibility is related 
to the transverse contributio»i2i2ii^. Specifically, the 
noninteracting Dirac model yields an orbital magnetic 
susceptibility Xmiq — > 0) oc (5(/i), in the long wave- 
length limit This is consistent with earlier results 
for graphite^'^, obtained using Wallace's two-dimensional 
band structure for a graphene layer-2d.. Such a finding 
would predict no response to a uniform, static mag- 
netic field, away from half-filling. This is of course par- 



tially compensated by a smearing of the (5-function at 
finite temperatures, already in the noninteracting limit. 
Still at zero temperature and in the noninteracting limit, 
one recovers a nonzero magnetic response also away 
from half-filling, when the honeycomb lattice structure 
is considered^. The effect of the interactions has been 
considered in Ref. [26l where it is shown that an inter- 
acting 2D Dirac electron liquid develops a magnetic re- 
sponse also at finite doping. Concerning the response to 
an external electromagnetic field, graphene is also unique 
among other conventional 2D electron systems, in that it 
has been predicted that it can sustain a transverse plas- 
mon mode^^, as a consequence of its double band struc- 
ture. More recently, such a transverse mode has been 
predicted also for bilayer graphene^l. 

Here, we will consider the effect of strain on the vari- 
ous electronic properties that may be described by linear 
response correlation functions. Indeed, a deformation of 
the lattice through the application of uniaxial strain or 
hydrostatic pressure is expected to produce modifications 
also in the electronic structure of graphene. Recently, it 
has been proposed that nanodevices based on graphene 
could be engineered on the basis of the expected strain- 
induced modifications of the deformed graphene sheet 
(origami electronics)SS,. This is made possible by the ex- 
ceptional mechanical properties of graphene, as is the 
case for other carbon compounds. For instance, despite 
its reduced dimensionality, graphene is characterized by 
a sizeable tensile strength and stiffness^^, with graphene 
sheets being capable to sustain elastic deformations as 
large as ~ 20%^"—. Larger strains would then induce 
a semimetal-to-semiconductor transition, with the open- 
ing of an energy gap^"— , and it has been demonstrated 
that such an effect critically depends on the direction of 
applied strain^i^. 

The paper is organized as follows. In Sec. |TT] we 
present our formalism for treating graphene under uni- 
axial strain, and derive our central result for a generic 
linear response function for a deformed graphene sheet. 
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This, in particular, applies to the density-density and 
current-current correlation functions. In Sec. lIIII we then 
study the electron polarization, with emphasis on the 
strain dependence of the plasmon dispersion relation and 
the conductivity. Specifically, we find a strain-induced 
anisotropic enhancement of the deviations from linearity 
of the transverse plasmonM, which should facilitate its 
experimental detection. In Sec. IIVI we derive the strain 
dependence of the magnetic and electric susceptibilities. 
Finally, we summarize our results in Sec. El 



II. MODEL 

The low-energy Hamiltonian for noninteracting quasi- 
particles around a Dirac point, say K, has the well-known 
linear form^ 

ij'"' = hvFO- ■ q, (1) 

where v-p is the Fermi velocity, <t = (ax,<Jy), with ax, 
ay Pauli matrices, and q is the wavevector displacement 
from the Dirac point one is referring to, i.e. k = K -I- q. 
Here and below, a superscript zero denotes absence of 
strain. The effect of strain is then that of modifying 

the lattice vectors as Si — (I + e) • S^'p {£ = 1,2,3), 
where 5^°^ = a(V3, l)/2, = a(-V3,l)/2, sf '^ = 
a(0,— 1) are the relaxed (unstrained) vectors connecting 
two nearest-neighbor (NN) carbon sites, with a = 1.42 A, 
the equilibrium C-C distance in a graphene sheet^, and 
e is the strain tensor— 

e=^e[{l-,y)l + {l + iy)Ai0)l (2) 

where 

A{9) = a,,e^'''''y = cos{29)a^ + sm{29)ax. (3) 

In Eq. ([2]), denotes the angle along which the strain is 
applied, with respect to the x axis in the lattice coordi- 
nate system, e is the strain modulus, and v is Poisson's 
ratio. While in the hydrostatic limit ly — —1 and e — el, 
in the case of graphene one has ly = 0.14, as determined 
from ab initio calculations^, to be compared with the 
known experimental value ly — 0.165 for graphite^. The 
special values 9 = and 9 = tt/6 refer to strain along the 
zig zag and armchair directions, respectively. 

The overall effect of a moderately low applied uniaxial 
strain on the low-energy Hamiltonian is that of shifting 
the location of the Dirac point in momentum space as 
K ^ and changing the shape of the Dirac cone, into 
a deformed one, with elliptical section'^. Such a picture 
applies for strain moduli below the value at which a gap 
opens in the energy spectrum, which takes place at e ~ 
20%^"—. In particular, setting k = k^ -I- q, with q 
measuring now the vector displacement from the shifted 
Dirac point, the Fermi velocity, defined as the slope of the 
Dirac cone in the direction of q, will now have anisotropic 



components C||1If, c±vf along the direction of applied 
strain and the direction orthogonal to it, respectively, 
with 

C|| = 1 - 2k£, (4a) 
c± = 1 + 2nue, (4b) 

where n = {a/2t)\dt/ da\ — ^ w 1.1 is related to the log- 
arithmic derivative of the nearest-neighbor hopping t at 
e = 0. Thus, the low-energy Hamiltonian around k^j 
maintains a linear form even in the presence of strain, 
and can still be written as 

H = hvFcr ■ q', (5) 

where now 

q' = R{9)S{e)R{^9)ci, (6) 

with R{9) the rotation matrix in the direction of applied 
strain, and S{e) = diag (c|| , c±) the matrix describing the 
deformation of the Dirac cone. Explicitly, for the com- 
pound transformation matrix R{9)S{£)R{~'9) mapping q 
onto q' one finds 

R{9)S{e)Ri-9) =I-2Ke. (7) 

A central result of the present work is that a similar 
correspondence holds between a generic linear response 
function x(q, w) under applied strain, with respect to its 
unstrained limit, x'^^H'^t'^)- This follows from the fact 
that any linear response function x(q, w) of a noninter- 
acting electron system can be expressed as an integral 
over the first Brillouin zone (IBZ) of a suitable matrix 
operator over pseudospins, which is itself a function of 
q. Such an operator then admits a unique expression in 
terms of the Pauli matrices ax, ay, az, and the identity 
matrix I = ctq. The simplest cases are then given by the 
density operator and the density current operator, which 
in reciprocal space read^^ 

p(°)(q) = E<-qJI*k, (8a) 

k 

(q) = -e^^F^I^k-qf^^^k, i^x,y, (8b) 

k 

respectively, where = (V'qA,V'qs)j ^ind V'qa destroys 
a quasiparticle with momentum q and pseudospin a = 
A,B, and summations run over the IBZ. While the den- 
sity operator does not change under applied strain, for 
the generic component of the density current operator 
one has 

J, = [I-2K£]„-jf\ (9) 

Here and below a summation will be understood over 
repeated indices (j = x,y). 

Defining now eigenvalues and eigenvectors in pseu- 
dospin space of the Hamiltonian with and without 
applied strain, Eqs. ^ and ©, as |q', A)(°) = 
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Mq' ^)^°^ and H\q, A) = £;Aq|q, A), respectively, with 
A a pseudospin index, it follows that both E\q and |q, A) 
under applied strain are mapped onto i?|ql and |q', A)^°\ 
respectively, where q' is given in terms of q by Eq. ©. 
Performing such a linear change of variables in the inte- 
gral defining the correlation function of interest, in the 
cases of the density-density and current-current correla- 
tion function, it follows therefore that 



n,,(q,c.) = [det5(£)]-^nW(q',c.), 
n,,(q,c^) = [det5(£)]-i 



(10a) 



X [I - 2KeUll'f^ iq', Lu) [I - 2Ke]fc/,10b) 



where detS'(£) = (1 - 2Ke){l + 2Kve). From Eq. ([TUa| . 
in the case of the density-density correlation function, it 
follows in particular that the effect of applied strain is 
that of transforming the momentum variable q into an 
'effective' one q', plus the introduction of an overall scale 
factor [det 5'(£)] , which is isotropic with respect with 
the strain direction. Such a scale factor is directly re- 
lated to the slope of the electronic density of states at 
the Fermi level. As is well known, this goes linearly with 
the chemical potential /i, and it has been shown that its 
steepness increases with increasing strain, for moderately 
low strain modulus^. In the case of the current-current 
correlation function, such an overall effect is then super- 
imposed to an anisotropic deformation, depending on the 
angle of applied strain, 0, as shown by Eq. (jlObp . 
Linearizing Eq. (|10p with respect to e, one finds 
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g4";(q,c.) 
dqh 



ShkQk, 



(11a) 



n,j{q,uj) = [l + 2K{l-u)e]ufj\q. 



dU^^{q,u^ 

—2k 



£hk<lk 



dqh 

2K{£,nW(q,c.)} 



(lib) 



where the curly brackets in the last term denote a matrix 
anticommutator. 



III. POLARIZATION 

A. Charge response: Plasmons and conductivity 

We now specifically turn to consider the density- 
density correlation function within linear response the- 
ory, i.e. the electron polarization npp(q, cj). Plasmon 
modes are then recovered as poles of the polarization, 
and the effect of strain on their dispersion relation has 
been studied in Refs. ITollTTI . In particular, by including 
local field effects, as a consequence of the two-band char- 
acter of the band structure of graphene, we found two 



main plasmon branches, the lower branch being char- 
acterized by the standard square-root dependence on q, 
at long wavelengths, as expected for a two-dimensional 
systemifl. 

In given limits, the asymptotic form of the noninteract- 
ing polarization in the absence of strain, say IIpp^ (q, uj) 
is known explicitly. For instance, in the long wavelength 
limit (q — )■ 0), one finds"^> 



nfj{q^0,uj) 



gs9vq^ 


'2fi 


h^log 


2n- 


Huj 




huj 


2Ai4 





-i — 0{huj 



2y) 



(12) 



where = gv ~ 2 take into account for spin and val- 
ley degeneracies, respectively. In other words, IIpp^ [q — > 
0,cli) = Z{uj)q'^, at a given uj, with the complex factor 
Z{uj) imphcitly defined by Eq. p^ . 

In the case of applied strain, but still in the nonin- 
teracting limit, this is then readily modified through the 
linearized Eq. (jlip . yielding 

npp{q-^0,uj) = [1 - 2k{1 + i^)e cos{20 ~ 2ip)] Z{uj)q^ , 

(13) 

where q = q{cosip,sm(p). Within the random phase ap- 
proximation (RPA), the interacting polarization reads 
npp(q,a;) = npp(q,w)/(l - V{q)Upp{q,uj)), where 
V{q) — e'^/{2ereoq) is the (bare) Coulombic electron- 
electron interaction, and is the dielectric constant of 
the medium. Solving for the plasmon dispersion relation, 
ReILpp{q,uj) = 0, at low energies one finds 



hujpi = \ - — /i [1 — k(1 + v)e cos(26' — 2</j)] 
V 27re 

= h£ji{ip)y/qa. 



(14) 



One thus finds that the prefactor uJi{ip) in the 
dependence is maximum [resp., minimum] for ip—6 ~ tt/2 
[ip—6 = 0], i.e. wavevector orthogonal [parallel] to the di- 
rection of applied strain. Correspondingly, one also finds 
for the imaginary part of the retarded polarizability along 
the low-energy plasmon branch 

Imnpp(q,w + iO+) = 



^ ^■^/x[l-K(l + z^)ecos(26l-2(/7)] 



X {qaf'/'^5{hjj - hjjpi{q)). 



(15) 



Therefore, one recovers a dependence of the plasmon 
spectral weight on the angle of applied strain, similar 
to that shown by Wi(iy9) in Eq. ([Ti| . 

Another quantity of interest which is related to the 
density-density correlation function is the optical con- 
ductivity, which can be obtained as 



T„„(a;) = — lim ^npp(q, w). 



(16) 
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Making use of Eq. and (|12p one therefore finds the 
optical conductivity in the presence of applied strain as 



aipip{uj) = (To [1 ~ 2k(1 + J^)£ cos(26' — 2(^)] 

(4 u i 
Q(hw - 2/i) + i--^ + - log 



2^ - 




) 


2Ai + 


hiL) 





,(17) 



where ctq = 7re^/2/i is proportional to the quantum of 
conductivity. In the hydrostatic limit, v — ~1, <t^^ does 
not depend on strain, as may be expected, as the un- 
strained relation does not contain the Fermi velocity. 

The above expression for the conductivity, Eq. (IT71) . 
can be exploited to study the strain dependence of the 
transverse electromagnetic mode, that has been recently 
predicted theoretically in graphene^, and in a graphene 
bilayer^-. In a 2D electron gas, the spectrum of electro- 
magnetic modes obeys the equations^ 
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FIG. 1: Showing deviations from linearity of the frequency of 
the transverse plasmon, Eq. (|19p . for q in the allowed range, 
for strain modulus e = 0.1, and strain direction ranging from 
Lp-e = Q (top) to V9 - 6^ = 7r/2 (bottom). 
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2eo C,(q,i^)c^ 



(18a) 
(18b) 



for the longitudinal and transverse plasmons, respec- 



tively, where C, {q,uj) 



(cj/c) , where c is the ve- 



locity of light in vacuum. While conventional 2D elec- 
tron systems cannot sustain a transverse electromagnetic 
mode, it has been predicted^"* that graphene can develop 
a transverse plasmon mode, as a consequence of a nega- 
tive imaginary part in the interband contribution to its 
optical conductivity, Eq. ([TT]). Its logarithmic divergence 
as ?za;//Lt — )■ 2 is in turn related to the discontinuous be- 
havior of the interband absorption of radiation at fre- 
quencies hw > 2/i. Such a feature is a generic conse- 
quence of causality, and is related through a Kramers- 
Kronig transformation to the step-like behaviour of the 
real part of the optical conductivity. This is in turn due 
to the existence of a Fermi surface, which is however ex- 
pected to be smeared at finite temperature, thus imply- 
ing the reduction of the logarithmic singularity into a 
pronounced (but finite) peak. 

Indeed, making use of Eq. (ITT)) in Eq. (llBap . one con- 
sistently recovers Eq. (|14l) for the longitudinal plasmons. 
On the other hand, substituting Eq. ([TTl) in Eq. (jl8bl) . 
one obtains the strain-dependence of the dispersion rela- 
tion of the transverse plasmon implicitly as 



he 



C{q,Lo) = {1 - 2k{1 + iy)e cos{29 - 2ip)) 





2fi + 




- 2 




2/i - 







(19) 



where a = e^/(47reofi.c) is the fine structure constant. 

Because of the small factor a in the left-hand side 
of Eq. p^ . the dispersion relation of such a transverse 
mode is close to the linear dispersion relation of the elec- 
tromagnetic radiation itself, uj — cq < 0. However, one 
may expect that applied strain enhances deviations from 
linearity (i.e., from the photon's dispersion relation), as a 



consequence of a strain-induced modification of the band 
dispersion. Fig. [1] shows indeed deviations from linearity, 
w — eg of the transverse plasmon, for q in the allowed 
range, for strain modulus e = 0.1, and strain direction 
< (fi — 9 < Tr/2. One finds indeed that, in the case of 
applied strain, deviations of the transverse plasmon dis- 
persion relation from that of the electromagnetic radia- 
tion become significant over a sufficiently wide window 
in hcq/fi < 2, especially when ip — = tt/2. 

Therefore, applied strain should help the experimental 
detection of this elusive collective mode. Indeed, the fact 
that the plasmon dispersion relation is close to the corre- 
sponding electromagnetic dispersion implies that such a 
transverse plasmon mode would have a marked photonic 
character, and a small plasmon linewidth would there- 
fore hinder its observation'*'^. On the other hand, at fi- 
nite temperature, the real part of the optical conductivity 
is nonzero also for hoj < 2/Lt, so that the transverse plas- 
mon mode does acquire a finite, albeit small, linewidtbi^. 
In particular, this applies to plasmon energies such that 
< 2fi — hu! < k^T. This is exactly where the plasmon 
dispersion relation deviates most from its photonic coun- 
terpart, the deviation being enhanced, and shifted away 
from the limiting case Huj = 2/i, in the case of applied 
strain, for q perpendicular to the strain direction. One 
therefore expects wavevectors of the order of Hcq < 2/i, 
or equivalently q/kp < 2vf/c <C 1, so that it is justified 
to employ Eq. 



B. Current response 

In the case of an applied vector field (e.g., an elec- 
tric field Ei), one may in general decompose the linear 
response function in a longitudinal and a transverse com- 
ponent as 



f \ I \ , fx 

X^J(q,a;) = — ^X||(9,w) + [Sij - —f 



X±{q,^), (20) 
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where q= |q|, for a homogeneous systemic. In particu- 
lar, in the case of the current-current correlation func- 
tion, the latter being proportional to the pseudospin- 
pseudospin counterpart, this can be further simplified as 

ng)(q,a.) = nf{q,u)5,j + hL"^ (g, c.)A,, (^), (21) 



Therefore, as a consequence of Eq. ([TO]), it follows that 
even though an unstrained system is characterized only 
by transverse response in the static limit--, it may de- 
velop a nonzero parallel response as a result of a strain- 
induced deformation. Making use of Eq. \llh\ , one finds 



where 



ni°)(q,^) = -[n}|"^(g,..)±nr(q,..)]. (22) 



(q, Lu) = Vif^ (q, w) - 2eK(l + v) HL''^ (g, w) cos(20 - 2^)% ^ nV^^ (g, w) Ay [B) + HL"^ (g, w) A, {ip + ^/4) sin(20 - 2(p) 



t(0) 



(0)/ 



t(0)/ 



-k[(1 - I/) + (1 + J/) cos(26' - 2(p)]e 



5n^^°^(q,tj) an^'^(g 



(0), 



9(7 



(5ij -f- q- 



(23) 



In the static limit (uj — 0), Eq. ([25)) can be further simplified, by considering the analytic result of Ref. with 



n[°^(q,0) = 0, and 



nf(q,o) 



gsgv& vy 
16hq 



1 - e {hvFq - 2^) ■ 



i^hv^q} hvpqV i^hvpq} 



- e (2/i - hvpq) 



(24) 



In particular, one recovers a vanishing response, 
nqq(q, 0) = 0, with Ilqq denoting the current-current 
correlation function for both vector potential and re- 
sponse field aligned with q, when q is aligned with the 
applied field also in the presence of strain, as expected in 
the static limit. 



IV. ELECTRIC AND MAGNETIC 
SUSCEPTIBILITIES 

A magnetic field applied in the direction perpendicular 
to the graphene plane can be described as 'Bcxtil) — 
Sea:t(q)z = «q X A, whcrc A = iiqy,-qx)Bcxt/q'^, in 
reciprocal space. The linear response to such a magnetic 
field is then given by a current Ji, which in turn produces 
a magnetization term = Xm'^cxt- 

In the case of a static, uniform magnetic field, oriented 
in the direction orthogonal to the graphene sheet, one is 
interested in the magnetic susceptibility defined as 

f dip 

XM = hni / — Xm(q,0). 



2tt 



(25) 



Making use of Eq. ([M[) . one obtains 

k(1 



XM 



lim 



r,2 



1 



i^)eq 



dq 



In the strained case, this reads 



XM = -Mo[l - 2k(1 - iy)e] 



6tt 



<5(m). 



(26) 



(27) 



One therefore obtains a qualitatively similar result to 
the case of undeformed graphene, treated within the 
Dirac approximation and neglecting the electron-electron 
interactioi*22i2^. On the other hand, applied strain causes 
a reduction of the magnetic response, Eq. ([2S[) . Al- 
though Eq. ([27p would imply no response to a static, 
uniform magnetic field away from half-filling, one ex- 
pects that finite-temperature effects would broaden the 
(5-function, already in the noninteracting limit. A qual- 
itatively similar smearing of the peak in the depen- 
dence on the chemical potential may also be induced by 
disorderiS. Still at zero temperature and in the nonin- 
teracting limit, one recovers a nonzero magnetic response 
also away from half-filling, when the honeycomb lattice 
structure is considered^^. The effect of the interactions 
has been considered in Ref. [2^ where it is shown that an 
interacting 2D Dirac electron liquid develops a magnetic 
response also at finite doping. 



An analogous procedure may be followed to derive the 
electric susceptibility Xe, entering the relationship i5E = 
Xe^cxt between the electric polarization and an external 
electric field. One is then interested in the static {lu = 0) 
limit of the density-density polarization. In the presence 
of applied strain, at arbitrary fi — hv-pk-p, using Eq. (|lip . 
one explicitly finds 
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9s9vq ^+ f ^IJ- 1 f.,f, o ^ 
— — Q{hvFq - 2fi) 



-4(1 - + (1 + ^) cos(2^? - 2^)]e|f^G< f 



Q{hvpq — 2/i), 



(28) 



wher( 



.17.18 



G^{x) = ±xVl - - arccosa;, |x| < 1. (29) 

In particular, at zero doping (p = 0, G^(0) = — 7r/2), 
one finds in general that 

Xe(q,0) = Viq)Upp{q,0). (30) 



It should be emphasized that, while Eq. (I30t describes 
the response of the system to a static electric field lying in 
the same graphene layer. More explicitly, in the undoped 
case, Eq. (15(11) reads 

Xe = limxe(q,0) = - ^^'^"^ 

X [1 + k(1 - v)£ - k(1 + j/)ecos(26i - 2(^)] ,(31) 

where Lp is the direction of the electric field on the 
graphene plane, thus showing that uniaxial strain in- 
troduces a modulation in the angle of applied strain. 
Moreover, comparison with Eq. ([26|) in the hydrostatic 
limit {v = —1) shows that strain enhances the electric 
response, while suppressing the magnetic one. 



V. CONCLUSIONS 

We have studied the dependence on applied uniax- 
ial strain of several linear-response electronic correla- 



tion functions of graphene. After deriving a general 
correspondence between strained and unstrained corre- 
lation functions, we have derived the strain dependence 
of the plasmon dispersion relation and of the optical 
conductivity. Specifically, we find that the prefactor in 
the y^-dependence of the plasmon frequency develops 
an anisotropic character, with maxima occurring when 
the wavevector is orthogonal to the direction of applied 
strain. Moreover, we derive a strain-induced anisotropic 
enhancement of the deviations from linearity of the re- 
cently predicted transverse plasmoni^, which should fa- 
cilitate its experimental detection in suitably strained 
graphene samples. Finally, we have compared and con- 
trasted the strain dependences of the magnetic and elec- 
tric susceptibilities, showing that strain enhances the re- 
sponse of strained graphene to an applied electric field, 
while suppressing the response to an magnetic field. 
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